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Introduction. The theoretical description of the ion transport in 
membrane systems in the galvanostatic mode is presented. A 
desalting channel of the electrodialysis apparatus is considered 
as a membrane system. The work objectives are the develop- 
ment and verification of a two-dimensional mathematical 
model of the stationary transport of salt ions in the desalting 
channel of the electrodialysis apparatus for the galvanostatic 
mode. 

Materials and Methods. A new model of ion transfer is pro- 
posed. It is based on the Nernst — Planck — Poisson equations 
for the electric potential and on the equation for the electric 
current stream function. A numerical solution to the boundary 
value model problem by the finite element method is obtained 
using the Comsol Multiphysics software package. 

Research Results. The developed mathematical model enables 
to describe the stationary transfer of binary salt ions in the 
desalting channel of the electrodialysis apparatus. Herewith, 
the violation of the solution electroneutrality and the formation 
of the dilated domain of space charge at overlimiting currents 
in the galvanostatic mode are considered. A good agreement 
between the physicochemical characteristics of the transfer 
calculated by the models for the galvanostatic and potenti- 
ostatic modes implies adequacy of the constructed model. 
Discussion and Conclusions. The developed model can inter- 
pret the experimental study results of ion transfer in membrane 
systems if this process takes place in the galvanostatic mode. 
Some electrokinetic processes are associated with the appear- 
ance of a dilated domain of space charge at overlimiting cur- 
rents. When describing the formation of this domain, it is pos- 
sible to find out how the processes dependent on it affect the 


ion transfer in the galvanostatic mode. 


Beedenue. CraTba ocBaljeHa TeopeTHuecKoMy OIHCaHHtoO 
Tipollecca MepeHoca HOHOB B MeMOpaHHBIX CHCTeMax B IaJib- 
BaHOCTaTHYeCcKOM pexuMe. B kayecTBe MeMOpaHHOM CHCTeMBI 
paccMaTpuBaeTca KaHall OOeccONMBAaHMA IIeKTPOWMaIM3HOTO 
ammapata. Llemu padorTsr: co3qaHve vu BepHduKallua WByMep- 
Hol MaTeMaTH4ecKol Moje cTal[MoHapHoro TepeHoca 
MOHOB COJIM B KaHasle oOeccONMBaHHA 9IeKTpOqMaM3Horo 
almapata Jd raibBaHOCTaTHYecKOro pexKuMa. 

Mamepuaaoi u memoovi. UpennoxweHa HOBad MOJeNb TepeHo- 
ca HOHOB. Ona ocHoOBaHa Ha CHcTeMe ypaBHeHuli Heppcta — 
IInaHka — Ilyaccona yy1a 9IeKTpHueckoro NOTeHIMala HU Ha 
ypaBHeHHu Id PyHKUWA 9eKTpHyeckoro ToKa. [lomyaeHo 
YucIeHHoe pellleHHe KpaecBol 3aya4uH MOJeIM MeTOOM KO- 
HeYHBIX 93JIEMCHTOB C TMOMOIIbIO MporpaMMHoro MakeTa 
Comsol Multiphysics. 

Pezyabmamoi uccaedoeanua. PaspaboTaHHad MaTeMaTH4eckad 
MOJeIb MO3BOJIAeT OMMCAaTb CTalMOHapHbI MepeHoc HOHOB 
OuHapHoli CoM B KaHasie oOeccOuMBaHHA 9IeKTpOTMasIM3HOro 
ammapata. [pu 9TOM y4HTbIBAaIOTCA HapyllieHHe IIeKTpoHel- 
TpauIbHOCTH pacTBopa H tbopMupoBaHHe paciiMpeHHoH oOsa- 
CTH TIpoOcTpaHcTBeHHOrO 3apsjla Ip CBepXMIpeeJIbHBIX TOKAX 
B TajIbBaHOCTaTHYecKOM pexuMe. OO alleKBaTHOCTH HOcTpo- 
eCHHOM MOjIeIM CBUeTeIbCTByeT XOpolliee coBlayqeHve Pu3u- 
KO-XHMHYeCCKHX XapaKTepHCTHK MepeHoca, paccuTaHHbIXx 10 
MOJICJIAM JIA PasIbBAaHOCTaTHYeCKOrO HM MOTeCHI[MOCTaTHueCKO- 
TO pexKHMOB. 

O6cyacdenue u 3akmouenun. Pa3spadoTaHHad MOJ[eJIb 1103B0- 


JIACT UHTCpIpeTHpOBaTh Pe3yJIbTaTbl SKCIICPHMCHTaJIBHbIX 
viccyIeqoBaHHi TlepeHoca HMOHOB B MeMOpaHHBIX CHCTeMax, 
ecIM ,aHHbIit Tipouwecc mWpoTekaeT B lajIbBaHOCTaTHYCCKOM 
pOoxkKUMe. Hekoroppie SIEKTPOKHHCTHYCCKHE TIPOWeCchl CBA3a- 
HbI C MOABJICHHCM paclimpeHHon oomacTH TIpOcTpaHCTBeHHOro 
3apaya TIpH CBepxiipeXesIBHBIX TOKAax. OnncpiBast cbopmuposa- 
Hue yka3aHHoli oosacTu, MO2XKHO BbIACHHTb, KAKHM o6pa30M 
3aBHCALIHe OT Hee TIpOWecchl BIIMAKOT Ha HWepeHOC HOHOB B 


TaJIBBaHOCTaTH4eCCKOM PpexKUME. 
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Introduction. Membrane systems form the basis for electrodialysis machines, nano- and microfluidic devices, 
which are used in water treatment, agricultural products (milk, wine, etc.) processing, performing chemical analysis, 
and in other areas of activity [1-4]. In numerous mathematical models of mass transfer processes in membrane systems 
for potentiostatic or potentiodynamic modes, the electrical mode is defined as a potential jump between two equipoten- 
tial planes parallel to the membranes. A detailed review of the papers devoted to modeling for the potentiostatic mode is 
presented in [5—7]. 

All the while, in the practice of electrodialysis, electrochemical characterization of membranes (chronopotenti- 
ometry, impedansometry, etc.), the galvanostatic mode is often used, in which a constant average current density is 
maintained at the interface. A huge amount of experimental data has been compiled about this mode. They must be in- 
terpreted [8—10]. Studies in the mathematical modeling of the galvanostatic mode are conducted in several directions. 

The first direction is the inverse problem method. As is clear from the name, this is about solving an inverse 
problem: for the specified current density at the “solution - membrane” interface, the corresponding potential jump is 
found, and then the problem for the potentiostatic mode is considered [11]. The low efficiency of this method is due to 
the fact that its implementation requires multiple solutions to the problem in potentiostatic mode for one given value of 
current density. 

The second direction is the decomposition method. In this case, the system of Nernst — Planck — Poisson equa- 
tions is replaced by a system of decomposition equations [12—16]. The assumption of a quasi-uniform charge distribu- 
tion enables to obtain a model for the galvanostatic mode in the approximation of Ohm's law [17-20]. 

The third approach can be called the direct method. In this case, the equation replacing the Poisson equation is 
derived for the current density in the desalination channel [21]. 

The galvanostatic mode can be described differently — by numerical solution to the Nernst — Planck — Poisson 
equations for an electric potential with a special boundary condition that allows the current density to be set as a param- 
eter specifying the electrical mode in the system. In [22, 23] for the one-dimensional case, the time derivative of the 
electric potential gradient was determined as an explicit function of the current density. This distinguishes the authors’ 
approaches from potentiostatic models in which the difference of potentials is set. 

This paper presents a stationary model of the ion transfer process in membrane systems for the galvanostatic 
mode. It is based on the Nernst-Planck-Poisson equation system with the boundary condition that enables to establish 
the current density as a parameter that sets the electric mode in the system. This solution is similar to [22, 23]. The dif- 
ference is that the proposed model is two-dimensional and considers the variability of the current density along the 
channel. 

Materials and Methods. The desalting channel of the electrodialysis apparatus (EDA) formed by the anion- 
exchange (AEM) and cation-exchange membranes (CEM), is taken to mean a membrane system. A binary electrolyte 
solution is pumped through it at V, average rate. 

In Fig. 1, x is a coordinate normal to the membrane surface varying from 0 (border with AEM) to h/ (border 


with CEM); y is a tangential coordinate to the surface of the membrane varying from 0 (channel entrance) to / (channel 
exit). 
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0 Entrance h «x 
Fig. 1. Scheme of EDA desalination channel. Concentration profiles of C,; cations (solid line) and C, anions (dot line), 


V forced flow rate are shown. 


(System of equations. Consider a two-dimensional stationary case of the system of equations describing the 
transfer of a binary electrolyte with no chemical reactions [24]: 


I; 2" pC vo=D VCXCs, i=1,2, (1) 
RT 

-divj,=0, i=1,2, (2) 

E€,Ad =—F(z,C, +2,C,), (3) 

12 i, 4+ 2545,)' (4) 


Here, i , D,,z; and C, are, respectively, flux, diffusion coefficient, charge number, and molar concentration of the i- 
th ion; is electric potential; V is flow rate of the electrolyte solution; ¢9 is electric constant; ¢, is relative dielectric 


constant of the electrolyte solution (assumed to be constant); / is current density; F is Faraday constant; R is gas con- 
stant; Tis absolute temperature; j,, j,,7, 6, C,,C, are unknown functions of x and y coordinates. 

The Nernst — Planck equations (1) describe the ion flux due to electromigration, diffusion, and convection; (2) is 
the material balance equation in the stationary case; (3) is Poisson’s equation for the electric field potential; (4) is cur- 
rent density in the electrolyte solution. We assume that the velocity distribution in the channel corresponds to the 
Poiseuille flow [24]: 





xX xX 
V, =0, 1, = 41 *). (5) 


To determine the unknown functions listed above, it is necessary to set boundary conditions for the system (1) 
- (3). Let us consider two electrical modes: potentiostatic, when the potential jump at the system boundaries is set con- 
stant, and galvanostatic, when the density of the current flowing through the interface is constant. 

Boundary conditions for modeling potentiostatic mode. We assume that the surfaces of ion-exchange mem- 
branes are equipotential. The system (1) - (4) includes the potential of the electric field only in the form of derivatives 


with respect to the spatial coordinate. In this case, only Ad = 6(h, y)— (0, y) potential jump is significant, where Ad is 
a known function, so we set, for example: 
(0,y)=0. (6) 
Then 
Ad = const . (7) 
The conditions (6) and (7) determine the potentiostatic mode. Other boundary conditions are given below. 


At the “AEM - solution” interface (x = 0), the concentration of coions (cations) is determined from the condi- 
tion of continuity of their flow at the “membrane — solution” boundary considering the selective AEM properties [16]:. 


(S Poe 2 (0,y)-“E ER (8) 


“ex. RT rT Ox z,FD, , 
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where 7; (i = 1, 2) are effective numbers of transfer of counterions in the membrane (CEM and AEM, respectively); 7; 
are numbers close to 1, and also for 7, = 1 ideally selective membrane, and the condition (8) turns into the condition of 


membrane impermeability for coions. 
The concentration of counterions (anions) depends on the exchange capacity of the CEM, which can be specified 


as: 
C, (0,y)=C,, =N,C). (9) 


a 


Here, N, constant shows how many times this concentration differs from the concentration in the volume of the solution 
[25]. 


At the “solution - CEM” (x = /) interface for ion concentrations, the conditions similar to those at the “AOM - solution” 
(x = 0) border, are accepted: 





C,(4,y)=C,, =NC,, (10) 
F 1-7, )L.(A, 
aC, 2.0, %\(1,y)=$ Ey) a1) 
ox RT Ox z,FD, 
At the channel entrance (v = 0) a uniform distribution of ion concentrations is assumed: 
C(e0) SC. 1212. (12) 


The condition for the potential of the electric field is obtained from the equations (1) and (4) considering the 
absence of current flow through /, input (x, 0, f) = 0: 











Ob(x, 0) . ual ai OC, (x, 0) +z,D, OC, (x, 0) (13) 
oy F(z, D, + 23D, )C, oy 7 oy 
At the channel exit (v = /), ions are freely carried by the solution flow: 
(i,j, )= (i,-——zD, C,Vo-D,VC,+CV)=CV,, i=1,2. (14) 


The system (14) also means that the sum of the diffusion and migration tangential components of the anions 


and cations flow is equal to 0: 
eg (Het =12. (15) 
Oy RT oy 
A “soft” condition is accepted for the potential, meaning that there are no sharp changes in the potential at the 
channel exit: 


i, 1) =0. (16) 
oy 


The boundary problem, which includes equations (1) - (4) and boundary conditions (6) - (16), simulates the po- 
tentiostatic mode, and its control key condition is (7). 

Boundary conditions for simulating the galvanostatic mode. Under the simulation of the galvanostatic mode, 
the condition (6) remains, and (7) should be replaced by the condition related to the specified value of i,, average cur- 


rent density at the “solution — CEM” interface (x =h ). 
To derive such a boundary condition, we substitute the relations (1) into (4) and express the gradient of the 
electric field potential: 
RT 


V 
$ F?(z/D,C,+2z5D,C,) 





(I +F(z,D,VC, +z,D,VC,)— F(z,C, + z,C,)V). (17) 





Assuming x =h in (17), we obtain the relation that connects the gradient of the electric potential with the given 
value of the current density at the boundary, that is, the boundary condition at the “solution - CEM” interface: 


I,+Fz,D, ae + Fz,D, a 





RT 
F? 2 DC, +25D,C, 


op 
Ox 











(Ay) 


In this case, 7, current density must satisfy the condition: 


lt ; 
: [LG =i. (19) 
0 
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The mathematical model of the galvanostatic mode consists of the system of equations (1) - (4). The boundary 
conditions (18) and (19) replace the condition (7). The other boundary conditions coincide with those for the potenti- 
ostatic mode. 

The potential jump in the galvanostatic mode is a calculated value. 

Transformation of the boundary conditions for simulating the galvanostatic mode. The condition (19) is incon- 
venient for a numerical solution, since it contains an integral. One of this condition conversion options is given below. 

In the stationary case, the current density is a solenoidal vector. Indeed, if we multiply (2) by z,and sum up, 


then div’ = 0. Consequently, there is such y function that 











lay, She, (20) 
ox? oy 
Using n function, the conditions (18) and (19) are rewritten as follows: 
P ea -S14 FD Ti+ FD, = 
p y x x 
h, hy), 21 
ae | ») F? z, D.C, +2;D,C, ( ») @1) 
lt 1; 1 ; 
—[1,(h,y) dy =—= | Sh, ») dy = (nh, - (0) =, (22) 
Ls 1+ oy l 
The equation (22) can be rewritten as: 
nh, !) _ nh, 0) = =i . (23) 


To close the expression system, it is necessary to obtain an equation for yn function. To this end, as in [15, 16], 


we inject a linear differential operator, which is a vortex function (rotor) in the two-dimensional case, for W arbitrary 


rW) = [= th) (24) 


two-dimensional vector: 





It is easy to check that: 
1) r(Vu) =0 for any uw smooth function; 


2) r(uW) =(Vu,W), +ur(W) for any u smooth function and any W smooth vector. 


Here, (Vu, W), = aw, - cw, is skew-symmetric scalar product of Vu and W vectors, moreover, (a,a), = Ois for 
x” y 
any da vector. 
Applying (24) to the equation for current density (4), we obtain: 
r(1) = Fzr(j,) + Fzr(ja). (25) 


Using the formula of flows (1), we obtain the ratio: 














2g F — 
rd, = “See ee) —Dr(V C,)+r(CV), i=12. (26) 
Hence, considering the properties of r operator: 
+ F = ae 
ri) =- pea, (VC,,Vo), +(VC,7) +CrV), i=1,2. (27) 
nk - Of, a ' 
Considering (27) and r(/) = s ae = A7, the equation (25) can be written as: 
Xx 
Age || eb a P|) a ay P| 
RT Ox ox ) oy oy Oy } Ox (28) 
OV, av. 
+F | z, a +2, ae Vi-F | z, eG +2, OG Vi+F (26; + B50) ) - ee : 
ox ox ) ° oy oy )~ ox Oy 
It follows from (28) and (20) that 1 function is determined up to a constant, therefore we can assume: 
n(h,0) = 0. (29) 


Then we obtain from (23) 
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nh, !) 7 it - (30) 
The conditions (29) and (30) are boundary for y function. 


Mathematical model for the galvanostatic mode in a dimensionless form. For a numerical study of boundary 
value problems, it is convenient to go to a dimensionless form. So we can simplify the equations and find out the actual 
number and set of parameters that determine the system behavior. Dimensionless variables describe a class of similar 
processes characterized by the same value of dimensionless numbers. 

Characteristic values describing the problem. When simulating mass transfer processes in the EDA desalina- 
tion chamber, a number of characteristic values are taken: 

- for spatial coordinates — h intermembrane distance; 


- for ion concentrations — C, volume concentration of the electrolyte; 
- for speed — V, average speed of the forced flow; 
- for diffusion coefficients D =D ,D , (z ,-z ,)/(D , z,-D , z ,) diffusion coefficient of the electrolyte; 
- for electric potential — ~, = RT/F thermal potential; 
- for current density — i, = FDC, /h value (analogue of the limiting density of diffusion current); 
- for ion flow — j, = DC, /h diffusion flow. 
Transition formulas. We translate the equations into a dimensionless form through the following relations ( (u) 


index denotes the dimensionless variants of quantities): 




















qi = am y= y y = a por V , CH= G; 21, 
h h h V, C, 
& 24 16 Ri ee ee D, ey 
o” = : T= : I. ni” = ’ fig =— Fh =1,2, pe papel 
by Io FDC, Jo D 
The system of equations in a dimensionless form has the form ((u) index is omitted for simplicity): 
j, =-z,D,C, Vo—D, VC, + PeCV, i=1,2, (32) 
-divj,=0, i=1,2, (33) 
eAd =—-(z,C, +2,C,), (34) 
An=- (20, Ie) = |e 7p Pi gep Oe) 
ox ox ) Oy oy Oy ) Ox (35) 
ov, OV. 
+Pe} z, as +Z, oe V,—Pe| z, CG +Z, ae V. + Pe(z,c, +2Z,C,) u ga Oe ; 
Ox ox ) ° oy oy ) ox Oy 
P=2j +2), (36) 


The system of equations (29) - (35) contains two dimensionless numbers: the Peclet number Pe=V,h/D and 
€=€.€,RT / (C,h°F’) . The physical meaning of ¢ parameter is that it is double square of dimensionless Debye length - 
2 
ft e=2(15/h) ISI: 
Estimation of the parameter values shows that under natural conditions for electrodialysis, the Peclet number 
has 10° —10° order, ¢ number has 10° — 10°’ order, that is, it can be considered a small parameter. 


For computational convenience, we transform the system of equations through plugging the flux density (32) 
in the equations (33) and (36): 


div(-z,D,C, Vp-D, VC, +PeCV)=0, i=1,2, (37) 


fF =)'z,(-z,D,C,Vo-D, VC, + PeCV). (38) 


2 
i 
=1 


Thus, the system of equations contains the following unknown x, y functions: C,,C,,o, /,,/,. The fields of 
C,, C, concentrations and @ potential are determined by solution to the equations (37), (34), respectively. /,, 7, current 


density components are calculated using (38). Speed distribution (5) in dimensionless form: 
V.=0, V, =6x(1-x). (39) 


Information technology, computer science, and management 


iN 
Ww 
— 


http://vestnik.donstu.ru 


aN 
(ee 


Vestnik of Don State Technical University. 2018. Vol. 18, no. 4, pp. 426-437. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk Joucxozo zocyOapcméennozo mexnuyeckozo ynuueepcumema. 2018. T. 18, Ne 4. C. 426-437. ISSN 1992-5980 eISSN 1992-6006 





The presence of ¢ small parameter in the Poisson equation (34) means that the boundary problem is singularly 
perturbed. This significantly complicates its numerical solution, since such problems are stiff [26]. @ potential of the 
electric field and C,, C, ion concentrations change very quickly in a narrow boundary layer whose thickness is equal to 
1, Debye length [5]. 

To solve this problem, it is advisable to compact the computational grid in the boundary layer and use special 


methods for solving stiff problems [26]. 
Boundary conditions in dimensionless form. At the “AEM — solution” interface (x = 0): 




















1-7, )I,(0, 
aa ogc (0,»)=! 2) FAC y) (40) 
ox ox z,D, 
C,(0,¥)=Nas (41) 
(0,¥) =0, (42) 
0 
“0, y)=0. (43) 
Ox 
At the “solution — CEM” interface (x = 1): 
C(Ly)=N, (44) 
aC, ab (1-7,)1,(Ly) 
+2,C, = |(1,y) =a 45 
ox es 2) ¥) z,D, ” 
ae -T14 2 Ft 2D, = 
yy x x 
—(1,y)=-| —— —————————_ |(L y), 46 
a ») z, DC, +25D,C, ( ») (6) 
on 
—(l,y)=0. 47 
Pa y) (47) 
At the channel entrance ( y = 0): 
C,(x,0)=1, i=1,2, (48) 
aD) ae: ia | easy, (49) 
oy z,D,+23;D, oy oy 
n(x,0) =0. (50) 
At the channel exit (y = /): 
eG, z,C, 2) s)=00=12 (51) 
oy oy 
(5 1) =0, (52) 
oy 
n(x,1)=—i,,1, i, = const . (53) 


After numerical calculation of the system (34), (35), (37) - (53), A@ potential jump in the desalination channel 


is determined by the formula: 
so 1 
Ag =—[o(Ly)ay. (54) 
0 


The numerical solution is found by the finite element method using the Comsol Multiphysics package on an 
uneven computational grid (the density of grid elements is increased at the “solution —- membrane” boundaries) [27]. 

Research Results. The calculations are performed fore =1,9-10~ , Pe = 2355, which corresponds to the fol- 
lowing values of the system parameters: 

- input concentration of the electrolyte solution of NaCl: Cy = 0.1 mol/m’; 

- temperature: T= 298 K; 

- diffusion coefficients of cations and anions, respectively: D, = 1.33-10° m’/s, D, = 2.05 - 10° m’/s; 

- numbers of counterions transfer in the membranes: T, = 0.972, T> = 1; 

- ion charge numbers: z; = 1, z. =—1; 

- ratio of the counterions concentration at the boundary with the membranes to its value at the channel en- 
trance: VN, =N,=1; 
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- channel width: h=10° m; 

- channel length: / = 2:10° m; 

- rate of the solution pumping: Vp =3,8-10° m/s. 

Fig. 2 shows the fields of C; u C2 concentrations, @ potential and y functions calculated at i4,=1,5ij current 
density, where i;;,, is the limiting current density determined by the Leveque formula (in dimensionless form) [28]: 


1 hy, ) 
i= 1,47) —2 0,2 |. (55) 
T,-t, ID 


1 





Here, t; = 0.395 is kation transport number in the solution [9]. 
i 


1| ™% 
































6) 





Fig. 2. Concentration fields of cations C, (a) and anions C2 (b), @ potential (c) and y function (d). Model calculation for 
galvanostatic mode at i,, =1.5i,,, current density. 
The forced flow acts in the channel, therefore in the areas near the membranes, the depletion of the ion concen- 
tration increases with distance from the channel entrance (along the direction tangential to the membrane surfaces). Ac- 
cordingly, the thickness of the expanded space-charge region increases along the channel (Fig. 3). 
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Fig. 3. Results of numerical calculation based on the model for galvanostatic (solid lines) and potentiostatic (dot lines) modes at 
i, =1,5i,;,, current density in cross sections y= 0,1/ (1), y=0,4/ (2), y=0,9/ (3): C; and C, concentration profiles (a); en- 


larged fragment of fig. 3, a (b); p = z,C, + z,C, space charge density (c) 
The volt-ampere characteristic (VAC, curve | in Fig. 4) is calculated on the basis of the galvanostatic model. 
In this case, the specified current density changed: i,, =0; 0,005i jin; ...3 1,5ijim. 


40 





30 
8 
20 
10 
0 10 20 30 40 
Ag 
Fig. 4. VACs calculated on the models for galvanostatic (1) and potentiostatic (2) modes; i}; limiting current density (3) determined 
by formula (55) 


The limiting current density is determined by the point of intersection of the tangents to the VAC in the initial 
part and on the inclined plateau [5]. Here it coincides with i;;,,, value estimated using the approximate Leveque formula 
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(55). Besides, the VAC was calculated on the basis of the potentiostatic model under a change in the potential jump Ao 
= 0; 0.4; ...; 40 (curve 2 in Fig. 4). As Fig. 4 shows, VAC | and 2 are congruent. Thus, there is a unique correspondence 
to each potential jump of a certain current density, and vice versa. 

Fig. 3 and 4 demonstrate a fairly good agreement of various physicochemical transport characteristics calculat- 
ed on the models for galvanostatic and potentiostatic modes. This proves the adequacy of the model of the transport in 
galvanostatic mode built by the authors. 

The computational costs of the proposed model were estimated as follows. The time spent on solving boundary 
value problems for the galvanostatic and potentiostatic modes was recorded. Herewith, the specified accuracy values of 
the calculations and the system parameters were the same, andi, = 1,5 7,,,.. Then time indicators were compared. Thus, 


it has been found that the calculation on the galvanostatic model requires 1.6 times more time. This is due to the fact 
that: 

- galvanostatic model contains an additional equation for determining the current density distribution; 

- for the potential at x = 1boundary, the condition of the second kind (46) is established. 

Complex potential of the electro-membrane system for the galvanostatic mode. In [25], a generalization of the 
impedance of an electrochemical system using 7 function introduced under the electroneutrality conditions was pro- 


posed. A similar generalization is acceptable in this case. Note that when using electrochemical impedance, an object is 
considered only as a “black box”, and its intrinsic properties are determined indirectly. The results obtained above ena- 
ble to introduce the concept of a complex potential of the electro-membrane system: P = $+i-n . Complex potential is 


a coordinate function of any point inside the object, therefore, unlike the electrochemical impedance, it allows us to 
explore the intrinsic properties of the object. 

Discussion and Conclusions. A new technique of mathematical simulation of the stationary process of mass 
transfer in the galvanostatic mode for membrane systems is described. In this case, the two-dimensional case is consid- 
ered with the use of a special boundary condition that enables to specify current density in the system. The equations for 
the current electric function are presented. The numerical solution results for potentiostatic and galvanostatic models are 
in good agreement. This shows the adequacy of the proposed transport model in the galvanostatic mode. 

The developed model makes it possible to interpret the experimental studies results of the ion transfer in mem- 
brane systems if this process proceeds in the galvanostatic mode. Some electrokinetic processes are associated with the 
appearance of a dilated space-charge region under over-limiting currents. The description of the formation of this region 
allows us to find out how the processes dependent on it affect the ion transfer in the galvanostatic mode. 
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